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Abstract 

The study of low density, ultracold atomic Fermi gases is a promis- 
ing avenue to understand fermion superfluidity from first principles. 
One technique currently used to bring Fermi gases in the degener- 
ate regime is sympathetic cooling through a reservoir made of an 
ultracold Bose gas. We discuss a proposal for trapping and cooling 
of two-species Fermi-Bose mixtures into optical dipole traps made 
from combinations of laser beams having two different wavelengths. 
In these bichromatic traps it is possible, by a proper choice of the 
relative laser powers, to selectively trap the two species in such a 
way that fermions experience a stronger confinement than bosons. 
As a consequence, a deep Fermi degeneracy can be reached having 
at the same time a softer degenerate regime for the Bose gas. This 
leads to an increase in the sympathetic cooling efficiency and allows 
for higher precision thermometry of the Fermi-Bose mixture. 
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1 Introduction 



Superfluidity and superconductivity are phenomena at the heart of quan- 
tum mechanics of many-body systems. Their importance is not limited 
to condensed matter physics and some of the concepts involved in their 
understanding have been seminal in other contexts, most notably in quan- 
tum field theory pQ. A wealth of experimental informations on macro- 
scopic quantum transport has been collected by studying the Bose and 
Fermi isotopes of helium in the superfluid phase and the electron liquids 
in superconductors With the advent of Bose-Einstein condensates of 
dilute gases consequent to the development of innovative cooling technolo- 
gies 0J |B] , the class of systems available for studying superfluidity has 
been enlarged to many other species, namely various isotopes and hyper- 
fine states of alkalis as well as atomic hydrogen, metastable helium, and 
ytterbium. Besides the richer spectrum of investigable atomic systems, in 
the case of dilute gases one can exploit their intrinsically slower dynam- 
ics - low densities imply weak average interactions - to study formation 
and decay of interesting structures like vortices, a possibility very hard to 
achieve in the denser liquid helium. Moreover, unlike liquid helium, the 
precision achievable with atomic physics experimental techniques, and the 
possibility to control the theoretical many-body approximations, are two 
other reasons which have contributed to the fast development of this sector 
at the borderline between atomic and condensed matter physics. 

After the pioneering observation of Bose-Einstein condensation (BEC) 
in dilute gases of 87 Rb 0, 23 Na [H], and 7 Li UJ, a lot of experimental 
and theoretical activities have been focused on the signatures of super- 
fluid behaviour in this novel low-density state of matter (see for instance 
^3^2^]). Examples are the formation of vortices by means of optome- 
chanical driving ^3] an d mechanical stirring |141 [T5] , the spectroscopy of 
scissor modes QB], the studies of superfluid flow and the related onset of 
a critical velocity ^| ^]. These experimental achievements have been 
complemented by relatively simple, first-principle theoretical studies with 
significant progress in the comprehension of longstanding issues in liquid 
4 He, like critical velocities and vortex formation ^5H2U]- Even more inter- 
esting is the study, still completely open from the experimental viewpoint, 
of superfluid features in Fermi dilute gases |21| . The presence of a super- 
fluid phase is expected to occur in the deep degenerate regime via a sort 
of atomic Cooper pairing |221 \'2'.'>\ . on the basis of qualitative analogies to 
the case of 3 He and, more in general, to high density electron liquids in 
superconductors. 

In this paper we discuss in some length a recent proposal to confine and 
cool two-species Fermi and Bose gases in a bichromatic optical dipole trap 
|24| . This configuration allows for selective trapping of the two species with 
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different trapping strengths. Since the confinement determines the degen- 
eracy conditions, a regime can be chosen such that the Fermi temperature 
is much larger than the Bose-Einstein condensation temperature. In this 
case, a deep Fermi degeneracy could be achieved before (or in proximity) of 
the BEC phase-transition for the Bose component, leading to various ad- 
vantages in the search for superfluidity in dilute fermions amidst a thermal 
(or thermally-dominated) Bose gas. This would provide an unprecedented 
situation as compared to the already available 3 Hc- 4 He mixture. In princi- 
ple, many Fermi-Bose mixtures can be studied. We will focus the discussion 
on the best two Bose coolers available, 23 Na and 87 Rb - the workhorses for 
Bose-Einstein condensation - and the only two stable Fermi isotopes for 
alkali atoms, 6 Li and 40 K. 

The paper is organized as follows. In Section 2 we briefly describe the 
experimental techniques specific to the trapping and cooling of Fermi dilute 
gases, and then give an updated overview of the current experimental ef- 
forts in reaching Fermi degeneracy. In Section 3 we introduce our proposal 
by discussing the conservative trapping features in various configurations 
and for diverse combinations of Fermi-Bose mixtures. In Section 4 we dis- 
cuss evaporative cooling for the boson species by giving a specific example 
of its dynamics in the case of one-color and two-color single-beam optical 
dipole traps. In Section 5 we deal with sympathetic cooling, with particular 
regard to the limitations to the minimum achievable temperatures induced 
by the heat capacities of the two species. In the same Section, we discuss 
qualitatively also some methods to evidence a possible superfluid state in 
the Fermi component, and comment on the advantages of having a ther- 
mal Bose cloud rather than the only Bose condensed component. General 
features of our proposal are summarized and discussed in the conclusions. 

2 Status of the experimental searches for su- 
perfluidity in dilute Fermi gases 

As a general consequence of the Heisenberg principle, quantum degeneracy 
occurs at temperatures very similar for bosons and fermions in presence 
(or absence) of external confining potentials with comparable strength. It 
is therefore natural to apply to the Fermi gases the successful trapping and 
cooling techniques already developed for bosons and culminating in the 
observation of a Bose condensed phase. However, there are also very strik- 
ing differences between bosons and fermions in the degenerate regime. For 
instance, fermions enter into a degenerate regime without the sharp phase 
transition characteristic of bosons. Also, when dealing with Fermi systems 
one has always to face the effects of the Pauli principle which freezes most 
of the available degrees of freedom. The effects of the Pauli principle are 
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particularly felt in all the current experimental efforts to achieve full de- 
generacy and to evidence a superfluid phase in Fermi systems confined by 
means of magnetic traps. With this trapping technique only spin-polarized 
Fermi gases can be confined and cooling is ultimately obtained - like in their 
bosonic counterparts - by the selective removal of the most energetic part 
of the atomic cloud allowing for rethermalization of the remaining fraction. 
Consequently, since the atoms are polarized in the same hyperfine state, 
the Pauli principle forbids the s-wave elastic scattering. Then rethermaliza- 
tion becomes inefficient at low temperature where the contribution of odd 
angular momentum states (like p-wave scattering) is strongly suppressed, 
and this limits the efficiency of direct evaporative cooling among fermions 
in the same hyperfine state. Two routes to overcome this limitation have 
been implemented, dual evaporative cooling of fermions prepared in two 
different hyperfine states, and sympathetic cooling with a Bose refrigerant. 
Even in these situations the Pauli principle gives limitations, as the elastic 
scattering between different hyperfine states is inhibited by the so-called 
Pauli blocking [25]: the atoms available to elastic scattering are limited to 
the Fermi surface and their number is directly proportional to T/Tp. The 
use of sympathetic cooling with a boson reservoir is instead obstacled by 
the superfluidity of the latter , the fact that the specific heat of the Bose 
gas quickly vanishes fot T < T c , and ultimately by Pauli blocking. As a 
matter of fact, the lowest temperatures presently achieved for fermions are 
limited in the 0.05-0.2 T F range |2Zl I2HI EH EDI EH E2I E3 EH ESI EH E3 

Concerning the trapping techniques, the use of magnetic traps gives 
limitations in the combinatorics of trappable hyperfine states and interferes 
with the use of tunable homogeneous magnetic fields required to enhance 
atomic scattering via Feshbach resonances as predicted in [35] (see also [3J5] 
for a recent review) and observed in various atomic systems |40l 1411 l4"2l I43| . 
Feshbach resonances provide a mechanism that could be crucial to identify 
signatures of superfluidity even at relatively large temperatures, the so- 
called resonant superfluidity [4*41 EE51 141)1 147] . 

Some of the above mentioned limitations can be overcome by using 
optical dipole traps @HJ EH1 E01 IS • Both different hyperfine states and 
arbitrary magnetic fields can be used in this case. Optical dipole traps have 
been pursued as a way to obtain quantum degeneracy with purely optical 
tools, avoiding the complications of magnetic trapping [52]. After studies 
on degenerate Bose gases generated in magnetic traps and then transferred 
into optical dipole traps [53] , both BEC [j2] E5] and Fermi degeneracy [2""j] 
have been achieved directly in all-optical traps. Also, preliminary studies 
of Feshbach resonances for fermions in an optical trap have been reported 
in |5S]. More recently, studies of strongly interacting Fermi gases have 
been reported in the case of Fermi-Bose mixtures [37] and two-component 
Fermi gases jS8J EH] . In the latter case the expansion dynamics has been 
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interpreted as a possible evidence of resonant superfluidity as predicted in 

The experimental situation, updated to the Summer 2003, is summa- 
rized in Table 1. The two stable fermionic species 6 Li and K have been 



Atomic species 


N F 


T/T F 


Ref. 


4U K(9/2,9/2) - 4U K(9/2,7/2) 


7 x 10 b 


0.50 


JILA1 [23 


6 Li(3/2,3/2) - 7 Li(2,2) 


1.4 x 10 5 


0.25 


Ricel[2I| 


6 Li(3/2,3/2) - 7 Li(2,2) 


4 x 10 3 


0.20 


ENS1 |2H1 


6 Li(l/2,l/2) - 6 Li(l/2,-l/2) 


10 5 


0.50 


Dukel P| 


6 Li(l/2,l/2) - 23 Na(l,-l) 


1.4 x 10 5 


0.50 


MITlgOI 


40 K(9/2,9/2) - 87 Rb(2,2) 


10 4 


0.30 


LENS [SI 


6 Li(l/2,l/2)- 6 Li(l/2,-l/2) 


1.6 x 10 5 


0.15 


Duke2 |S21 


40 K(9/2, -9/2) - 40 K(9/2, -5/2) 


1.1 x 10 6 


0.21 


JILA2 [33 


6 Li(3/2,3/2) - 7 Li(2,2) 


7 x 10 7 


0.10 


Rice2|331 


6 Li(l/2,l/2) - 23 Na(2,2) 


7 x 10 7 


0.05 


MIT2|3SI 


6 Li(l/2,-l/2)- 6 Li(l/2,l/2) 


8 x 10 4 


0.43 


ENS2~ 36 



Table 1: Status of the experimental studies of degenerate Fermi dilute 
gases. For each laboratory the various trapped species are reported speci- 
fying the particular hyperfine state (-F, tuf), the number of fermions at the 
final stage of cooling N Fl the lowest reached temperature ratio T/Tp, and 
the related reference. The quoted experiments at JILA and Duke University 
make use of dual evaporative cooling of two hyperfine states of magnetically 
trapped potassium and optically trapped lithium, respectively. The other 
experiments exploit sympathetic cooling with bosonic reservoirs, Li, 23 Na, 
and 87 Rb. More recently, sympathetic cooling of 40 K with 87 Rb has been 
also pursued at JILA |37j. In subsequent work the Duke University group 
has reached deeper degeneracy and observed an anisotropic expansion of 
the Fermi gas, which could be interpreted as an evidence of superfluid be- 
havior of the Fermi gas [581 159) . More recent experiments use Feshbach 
resonances, resulting in high-efficiency production of potassium or lithium 
ultracold molecules. In the experiment described in [32] the Fermi-Bosc 
mixture is an intermediate stage, then all the bosons are removed and the 
fermions are prepared in an incoherent mixture of equal populations in the 
(1/2,-1/2) and (1/2,1/2) states. 

cooled by using evaporative cooling between two hyperfine states (JILA 
and Duke University) or through sympathetic cooling with Bose conden- 
sates of 7 Li, 23 Na, and 87 Rb (Rice-ENS, MIT, and LENS respectively). It is 
evident that, despite of very different trapping and cooling techniques, the 
lowest degeneracy parameter T/Tf obtained by using sympathetic cooling 
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is around 0.2 (for the particular case of MIT2 see a detailed discussion in 
Section 5). The existence of this sort of lower bound for T/Tp can be un- 
derstood semiquantitatively in the following way |SJ . In the latest stage of 
evaporative cooling the number of Bose atoms becomes of the same order 
of magnitude of the number of Fermi atoms (in some hyperfine states of 
7 Li, due to the negative scattering length, there is also a theoretical upper 
limit to the number of atoms in the condensed phase, see On the other 
hand the boson specific heat scales as T 3 below the BEC transition, while 
the Fermi atoms have specific heat scaling as T. The specific heat curves 
of bosons and fermions, assuming for simplicity equal masses for the two 
species, intersect each other at a temperature T* ~ 0.5 T c ~ 0.25 Tp. Below 
T* sympathetic cooling becomes very inefficient. This simple explanation 
gives also an hint on how to overcome the limitation. If one could engineer 
the trapping potential in such a way that T c -C Tp, the Bose gas would 
preserve enough thermal capacity to drive the sympathetic cooling even in 
a deep degenerate regime for the Fermi component. A first attempt in this 
direction can be found in |62j . where adiabatic compression was proposed 
in an optical dipole trap superimposed to an already confining magnetic 
trap. This should allow a small fraction of the Fermi atoms to experience 
a very tight confinement potential, thus enhancing the Fermi temperature. 
More recently, we have proposed an alternative solution in the context of 
pure optical trapping [23], by using a two-color optical dipole trap which 
confines both Fermi and Bose gases with different strengths. It is the pur- 
pose of the next Section to discuss in detail the static confinement features 
of this class of atomic traps. 

3 Two-color optical dipole traps 

Let us start our analysis from the conditions of degeneracy for Fermi and 
Bose gases confined in harmonic traps. The Fermi and Bose temperatures 
for dilute atomic clouds trapped by harmonic potentials can be written as: 

T F = 6 1/3 HuJiN^k^ 1 ~ 1.82 hujiN^k^ 1 (1) 

T c = C(3)- 1/3 tkJ h N^ 3 k^ ~ 0.94 fiu h Nl /3 k^ (2) 

with ujf = (ujfxUjfyLUfz) 1 / 3 and lu^ = (LU\ >x uj\>yUj\> z ) 1 t' i being the geometri- 
cal average of the angular trapping frequencies in the three directions for 
fermions and bosons, N{ and N\> the number of atoms of the Fermi and Bose 
gases, and h and k-Q the Planck and Boltzmann constants, respectively. 

Besides the small difference in the prefactor, the degenerate tempera- 
tures for Fermi and Bose atoms per unit of atom are similar in traps with 
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Figure 1: Combined optical dipole trap for two-species mixtures. Case of 
a single beam configuration (left) with the two color beams propagating in 
a coaxial fashion, and of a coplanar crossed-beam dipole trap (right) with 
relative angle between the two beam pairs 6 = 0. The lenses are assumed 
to be achromatic to ensure a common focus for the collimated beams at 
different wavelengths. 



the same angular trapping frequencies for the two species. This is indeed 
the case of magnetic traps: since the magnetic moments of the alkali- 
metals are very similar, the only difference in the trapping strengths is due 
to their different masses, with the angular trapping frequencies scaling as 
ujf/ujb ~ (rob/mf) 1 / 2 . The situation may change in optical dipole traps 
where the confinement is dictated by the detunings of the laser beams with 
respect to the atomic transitions and by the beam intensities. Situations 
for which Fermi degeneracy is reached before BEC (i.e. T-p > T c ) are 
therefore viable, provided that fermions and bosons have different atomic 
transitions. This will restrict our analysis to two-species mixtures, since 
the isotopic shifts are usually not enough to ensure selective trapping in 
single species Fermi-Bose mixtures without incurring in prohibitive heating 
due to residual Rayleigh scattering. 

In discussing optical dipole traps, one can consider either single beam 
or crossed-beam configurations [52]) see Figure 1. The former has the 
advantage of being simpler with fewer experimental problems of loading 
and alignment, the latter gives rise to a more isotropic confinement. In the 
following, we will discuss in detail the crossed-beam geometry. Results for 
the single beam configuration will be obtained as a particular case in which 
one of the two beams is turned off. 

In the crossed-beam configuration, a pair of laser beams red-detuned 
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with respect to the atomic transitions and focused on the center of the pre- 
existing trapping potential (for instance the one generated by the magneto- 
optical trap typically used for precooling the atomic clouds), gives an ef- 
fective attractive potential for both the species. This attractive potential 
is partially balanced by a second pair of mutually orthogonal blue-detuned 
laser beams, acting for instance along the same plane formed by the red- 
detuned beams and forming with the latter an angle 9. This second pair 
of beams gives rise to a repulsive potential and, by a proper choice of its 
detuning and power, provides a selective deconfinement for the two species. 

The potential generated by the dynamical Stark effect, felt by an atom 
of species a [a — b for bosons, a = f for fcrmions) whose atomic transition 
wavelength and linewidth are respectively A Q and r Q , and due to the laser 
beams i of wavelength Ai and intensity Ii (i = 1, 2 for the red-detuned and 
blue-detuned laser beams, respectively), is 0Hj 



where f2 Q = 2nc/X a , ^; = 2irc/Ai, and I^, at is the saturation intensity 
for the atomic transition, expressed in terms of the former quantities as 
jsat _ fi£l^Y a /12irc 2 . It is easy to recognize that the potential energy 
has the same sign of the laser intensity if fli > f2 a , i.e. for blue-detuned 
light. Then the atoms, trying to minimize their potential energy, move 
towards the regions of space with minimum light intensity, and therefore 
are expelled by the laser beam. The opposite occurs for red-detuned light, 
for which the atoms are attracted in the regions of maximum light intensity. 

The laser intensity is the incoherent sum of the intensities of each pair of 
beams (such incoherent sum can be obtained by orthogonal polarizations or 
by slight detuning of the beams within each pair) , all focused at the origin 
of the trap. If we assume that the red-detuned beams propagate along 
the x-y axes while the blue-detuned beams propagate along the axes 
possibly rotated by an angle 9, i.e. £ = x cos 9 + y sin 9, r] = y cos 9 — x sin 9, 
with < 9 < 7r/4, the intensities can be written as 



where Pi is the laser power, Wi is the 1/e 2 beam waist radius, and R{ = 
muf/Xi the Rayleigh range. The quantity 2Pi/Trwj = 7^(0,0,0) represents 




(3) 
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the peak laser intensity due to each pair of beams, obtained in the focal 
point. 

The overall potential felt by the fermions (bosons) is Uf = U[ + U\ 
(Ub = Ui + t/k)- For a proper choice of the laser powers, these potentials 
are energy wells with depths AUf and AUb which constitute the confining 
energies of the fermionic and bosonic species. The curvatures around the 
minimum of these potentials determine the trapping frequencies LOf and u>b 
and, consequently, the critical temperatures Tp and T c . Up to quadratic 
terms, the potentials Uf and Ub are invariant under rotations around the 
minimum so that the trapping frequencies do not depend on the angle 9. 
By neglecting the terms {Xij-nwi) 2 « 1, we find 



Uaz / h fkfP 1 k%P 2 

u ax = u av = = y — + ] , 

where the constants kf have been introduced as 

1 1 



(7) 



Analogous expressions can be obtained in the case of a single-beam ge- 
ometry. For instance, by considering a single beam propagating along the 
x-axis we have 

/ h fkfP, k°K\ 
ax ~ V vm a \wlRl + w\R\ ) ' K ) 



h fkfPf k%P 2 



Jql V 



7Tm a \ W? Wn 



(9) 



Since typically Wi » A^ and, therefore, Ri ^> Wi, the confinement along 
the x-direction is weakened leading, with respect to the crossed-beam con- 
figuration, to an overall reduction of the average angular frequencies for 
the trapped species. In the particular case of equal ratios between wave- 
lengths and waists for the two colors, \i/wi — X2/W2 = X/w, this reduction 
amounts to (X/2ttw) 1 ^ 3 . 

In Figure 2 we show the dependence of the optical potential upon the 
radial, x or y, and axial, z, directions in the two-beam geometry for the case 
of the 6 Li- 23 Na mixture. When the ratio between the power of the blue- 
and red-detuned lasers is small (P2/P1 = 0.05, upper panels) the deforma- 
tions induced by the blue-detuned beam are negligible. The difference in 
curvature between the two species is mainly attributable to the difference 
in mass m a and detuning with respect to the laser beam, both playing a 
role in favoring T c smaller than Tp. Also, the depth of the potential well for 
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Figure 2: Potential energies in the case of the 6 Li- 23 Na mixture in a bichro- 
matic optical dipole trap in a crossed-beam geometry formed by focusing 
a Nd:YAG laser beam (Ai =1064 nm) and a blue-detuned laser beam at 
its second harmonic (A2 =532 nm). In the two upper panels the potential 
energy for the atoms of the two species is shown along the radial, x or y, 
and axial, z, directions for a power ratio P2/P1 = 0.05. The four curves 
in the radial directions refer to various angles (9 = 0, tv/16, 7r/8, and 7r/4 
from top to bottom) between the two pairs of beams. In the two lower 
panels the case of a larger power ratio, P2/P1 = 0.25, is considered. The 
beam waists are w\ — w% — 8 /im, and P\ = 1 W. Note that while the 
potential vanishes with a purely Gaussian behaviour in the axial direction 
z, along x (and y) first it varies in a Gaussian way (to half of its peak 
value for = 0) and finally vanishes as a Lorentzian (not shown). The 
Gaussian and Lorentzian widths are determined by the beam waists and 
the Rayleigh ranges, respectively. 
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mixture A f (nm) T f (MHz) A b (ran) T h (MHz) A 2 (nm) 



e Li-^Na 


671 


2tt x 5.9 


589 


2tt x 9.8 


532 


6 Li- 87 Rb 


671 


2tt x 5.9 


780 


2tt x 5.98 


740 


40 K- 23 Na 


767 


2tt x 6.09 


589 


2tt x 9.8 


532 


40 K- 87 Rb 


767 


2tt x 6.09 


780 


2tt x 5.98 


773.5 



Table 2: Wavelenghts and linewidths of the atomic transitions for the 
fermion-boson mixtures considered in the text and wavelength of the corre- 
sponding deconfining laser with a Nd:YAG laser at Ai = 1064 nm producing 
the primary trapping potential. 

the bosonic species is smaller than that for the fermions. This makes possi- 
ble to exploit evaporative cooling without appreciable interference from the 
Fermi cloud. It is also evident that the stronger confinement in the radial 
direction is achieved for coaxial beams, i.e. 9 = 0. The axial confinement 
along the z axis is instead unaffected by the rotation angle between the 
beam pairs. In the case of a strong perturbation, lower panels of Fig. 2 
where P2/P1 = 0.25, the bosonic species is much less confined and there is 
a strong difference in the curvature of the energy potential with respect to 
the fermionic species. 

In both crossed- or single-beam configurations, the confinement energy 
of the boson species vanishes as the power ratio P2/P1 approaches a critical 
value. Approximately, this can be obtained by Eq. 10 or Eqs. I|8I9|) as the 
P2/P1 value at which the average trapping frequency of bosons vanishes. 
In both cases, we find 



The behavior of the trapping frequency ratio uJi/ujh as well as of the 
trapping energies A?7f and AUb as a function of the power ratio P2/P1 is 
shown in Figs. 3 and 4 also for different Fermi-Bose mixtures, obtained by 
using fermionic 40 K and bosonic 87 Rb, as summarized in Table 2. The 
trapping frequency ratio for the single-beam geometry is always reduced 
with respect to the corresponding crossed-beam case at the same power 
ratio P2/P1. This decrease, which is a consequence of the weaker confine- 
ment along the direction of the laser beam, can be compensated by a larger 
power ratio, provided that the beam intensities are adequately stabilized. 
From the upper panels in Figure 4 we see that, with the chosen values of 
the blue-detuned laser wavelength, the bosonic sodium atoms are always 
less confined than the fermionic atoms. This allows one to apply evapora- 
tive cooling techniques for the bosonic species while having negligible losses 
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Figure 3: Selective trapping of bosons and fermions. The ratio between 
the average trapping frequencies for the fermionic and bosonic species is 
shown versus the beam power ratio between the blue- and red-detuned 
lasers. The cases of single beam (dashed lines) and crossed beam (solid 
lines) geometries are depicted for the Li-Na, Li-Rb, K-Na, and K-Rb mix- 
tures. The vertical lines are the critical values P2/P1 given by Eq. (JTUJ. 
Parameters chosen according to the values in Table 2. 
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Figure 4: Selective trapping of bosons and fermions. Confining energy per 
unit of infrared laser power as a function of the beam power ratio between 
the blue- and red-detuned lasers. The confining energy is evaluated with 
respect to the region of space delimited by the waist size (see also note 
in the caption of Figure 2). Continuous lines refer to the crossed-beam 
configuration, while dashed lines are for the single-beam optical dipole 
trap. 
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in the number of fermions. For the mixtures utilizing 87 Rb, the confining 
energy is initially higher for the Bose species (especially for the Li-Rb mix- 
tures) , however the fermionic confining energy increases as the laser power 
P 2 is increased, since the potential created by the blue-detuned light is in 
this case attractive for fermions (A2 > Af). For these mixtures, therefore, 
the addition of the laser light deconfining the bosons improves at the same 
time the confinement features of the fermionic species and minimizes the 
losses of Fermi atoms at the beginning of evaporation. Note also that for 
all the mixtures considered, within a very good approximation the bosonic 
confinement energies AC/b in Figure 4 vanish at the same critical values 
P2/P1 at which the ratios u>{/u)b in Figure 3 diverge. 

Some further comments are in order. In the case of equal waists, the 
potential felt by the bosons is exactly zero at the critical power ratio and 
changes curvature at larger values giving rise to deconfinemcnt, see mid- 
dle panel of Figure 5. By choosing W2 < Wi (top panel of Figure 5), 
the potential becomes bistable for P2/P1 around the critical value. Then, 
the blue-detuned laser can be used to produce tailored bistable potentials, 
providing an alternative route to the recently demonstrated all-magnetic 
bistable potentials |63|. For w 2 > ^i (bottom panel of Figure 5), the poten- 
tial is always globally deconfining but a local minimum around the origin 
is assured if P2/P1 is close to the critical value. The strong dependence of 
the critical power ratio upon the beam waists can be used to reduce the 
amount of blue-detuned light necessary to approach the targeted values of 
ujf/ajb- However, based on the above mentioned considerations, there is a 
tradeoff since we may have also a change of the potential shape. 

A possible technical issue is the demand for an initial large laser power 
especially since it is difficult to get sub-Doppler cooling for Li and K due 
to the small hyperfine splittings of the p 3 / 2 atomic level jBJ (see however 
for a successful demonstration of sub-Doppler cooling of 40 K). Up to 
now, Nd:YAG laser powers with built-in crystals for frequency-doubling are 
limited to about 1 W. This corresponds to order of 500 /iK for the initial 
potential depth, i.e. four times the initial temperature of the cloud if the 
latter is transferred from a Doppler-limited magneto-optical trap. As we 
will discuss in the next Section, this could be an issue for starting efficient 
evaporative cooling. Besides awaiting progress in the power deliverable by 
semiconductor-based lasers, one can use independent systems for red- and 
blue-detuned light (this seems anyway unavoidable for mixtures utilizing 
87 Rb as the Bose cooler). High power, far detuned CO2 lasers are also avail- 
able to produce quasi-static optical dipole traps H>3 . This technique 
has been recently demonstrated for cooling fermionic lithium at degener- 
acy temperature |29j after efficient loading from a magneto-optical trap 
jHHIEni- C0 2 lasers allow for large powers (order of 100-200 W) and, since 
they are well detuned from the atomic resonances, for very small heating 
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Figure 5: Bosonic potential energy along the z-axis in the case of the 
6 Li- 23 Na mixture in a crossed-beam bichromatic optical dipole trap for 
PijP\ 5% smaller (solid), equal (dashed), and 5% larger (dot-dashed) 
than P-2,jP\\ C nt- The three panels refer to waist ratios wijw\ — 0.9 (top), 
wijw\ = 1 (middle), and wijw\ = 1.1 (bottom). 
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Figure 6: Selective trapping with a CO2 laser in the case of the K-Na 
mixture. Trapping frequency ratio (left) and confinement energy (right) 
versus the parameter P-i/Pi. The dashed curves refer to the single-beam 
configuration, the continuous ones to the crossed-beam geometry. The 
different critical values of P2 /Pi obtained from the left and right panels are 
due to the fact that the harmonic approximation used to evaluate uj{/lo^ 
fails for P2/P1 large enough (the potential becomes bistable since u>2 < w\). 
The waist of the CO2 laser beam is assumed to be w\ = 50 //m, while 
w 2 = 10 fim. 



due to photon scattering |7D . In Figure 6 we show the selective trapping 
features of a K-Na mixture with the Nd:YAG laser replaced by a CO2 laser 
emitting at 10.6 /xm. At such large wavelengths the optical potential is 
well approximated by the induced dipole interaction U = —a g E 2 /2, where 
a g is the ground state polarizability and E 2 is the time-averaged squared 
electric field In terms of the laser intensity the optical potential is 

U{x,y,z) = I(x,y,z) (11) 

c 

Due to the high static atomic polarizability of rubidium (ofg b = 47.3 x 
10~ 24 cm 3 ), this bosonic species is more strongly confined than fermionic 
potassium (af = 43.4 x 10~ 24 cm 3 ) and lithium (a^ = 24.3 x 10~ 24 cm 3 ), 
therefore ruling out evaporation as a possible cooling technique for K-Rb 
and Li-Rb mixtures. Also, the combination Li-Na has to be discarded due 
to the small static atomic polarizability of sodium (a g a = 24.08 x 10~ 24 
cm 3 ). Thus, a CO2 laser can be effectively used for a red-detuned attractive 
potential only in the case of the 40 K- 23 Na mixture considered in Figure 6. 
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4 Evaporative cooling 



So far we have considered the static trapping features of the two-color 
optical dipole configuration schematized through a time-independent po- 
tential. However, to reach quantum degeneracy, the phase space density 
has to be increased, e.g. by cooling down the atomic sample. In the config- 
uration proposed here this is obtained through two continuous processes: 
forced and selective removal of bosonic atoms with re-thermalization of the 
surviving component (evaporative cooling) [6], and thermalization of the 
fermionic species to the temperature of the Bose gas (sympathetic cooling) . 
These cooling processes have various limitations, the most obvious one be- 
ing the reduced heat capacity of the bosonic sample undergoing a contin- 
uous decrease of atoms. Moreover, concurrent heating sources will limit 
the ultimate reachable temperature, most notably the residual Rayleigh 
scattering from the trapping beams. In this Section we discuss to some 
extent these issues and, although not pretending to go into full details on 
the many possible experimental schemes, we give estimates of the relevant 
parameters involved in the cooling dynamics for a specific example. 

Evaporative cooling in an optical dipole trap has been the subject of 
various experimental and theoretical studies |521 1541 ITT] . The atomic evap- 
oration rate in finite-depth traps is exponentially dependent upon the ratio 
between the potential depth of the trap and the average thermal energy 
of the atomic cloud, AU/IcbT, since the number of atoms in the tail of 
the Boltzmann distribution scales as exp(— AU/ksT). Forced evaporation 
is, therefore, necessary to maintain a significant evaporation rate, which 
should be otherwise exponentially quenched by reducing the cloud temper- 
ature. It can be demonstrated |71j that all the relevant quantities involved 
in forced evaporative cooling in an optical dipole trap scale with some 
power of the confining energy, provided that the ratio between the latter 
and the cloud temperature is constant in time AU /k^T = 77. In this case, 
the confining energy has a time dependence of the form 



where AU[ is the initial potential depth, ejj — — 2(f/ — 3)/f]', and r" 1 = 
(2/3)7/(77 — 4) exp(— 77)71, with -q' = 77 + (77 — 5)/(?7 — 4) and 7; being the 
initial value of the elastic collision rate 7. Similar laws hold for the time 
dependence of the number of particles N, temperature T, phase-space den- 
sity p, and elastic collision rate 7 with corresponding exponents ejv, £t, £p, 
and e 7 . 

With respect to the simplest case of a single Bose species undergoing 
evaporative cooling in an optical dipole trap we have two differences. First, 
there is a mixture also containing the Fermi atoms. Their presence does 




(12) 
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not significantly affect evaporative cooling of the Bose gas if the potential 
energy depth ALVf is much larger than AUb, a condition well satisfied in the 
situations we consider - see Figure 4. Also, in order for the Bose gas to act 
as a cooler its heat capacity must be much larger than the heat capacity 
of the fermionic species. This is possible, at least in the earlier stages of 
evaporation before entering the degenerate regime, by assuming a number 
of bosons much larger than the number of fermions. In bichromatic traps, 
a further problem may arise due to the presence of the blue-detuned laser 
which weaken the confinement of both species, particularly the bosonic one. 
As the atomic densities are decreased, all the density-dependent scattering 
rates, in particular the elastic one which is crucial for thermalization of the 
surviving atoms, are also suppressed [72] . 

The lower densities and the slower thermalization at lower temperatures 
makes more difficult to achieve high phase-space densities, an issue which 
has limited for various years the effort to make an all-optical formation of 
Bose-Einstein condensates [73]. From this point of view, in a bichromatic 
trap it is convenient to turn on the blue-detuned beam only when it can 
give an effective advantage, i.e. close to the Fermi degenerate regime. 
Once on, the power of the blue-detuned beam must be chosen in such a 
way to maintain a fixed ratio with respect to that of the red-detuned beam. 
Meanwhile the power of the red-detuned beam is decreased continuously 
to allow efficient forced evaporation of the bosonic species. An example of 
this cooling strategy is shown in Figure 7. 

The constant power ratio PijP\ must be chosen carefully as a com- 
promise between increasing the fermion-to-boson trapping frequency ratio 
and not decreasing too much the absolute frequencies, as this will decrease 
all the elastic scattering rates crucial for interspecies and intraspecies ther- 
malization. In fact, the elastic collision rate for indistinguishable bosons is 
given by 

7 = N h m h Lulcj h /2Tr 2 k B T, (13) 

where a\, — 8ira^ is the boson elastic cross-section expressed in terms of 
the boson s-wave elastic scattering length a^. The value of 7 during the 
cooling strategies described here is shown in the bottom panels of Figure 
7. 

If the technical noise and heating sources are properly reduced and the 
residual background pressure is low enough, the ultimate heating source 
is set by residual Rayleigh scattering from the laser beams. This depends 
on various parameters, most notably the laser intensity and the detuning 
of the laser beam from the atomic resonance, and the corresponding rate 
reads 

i?(zvz) = *$l V (^— + _J_V J (*,y,*) (14) 
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Figure 7: Evaporative cooling strategies in Dichromatic optical dipole traps. 
The time dependence of the laser powers (top panels), the averaged fre- 
quencies (center panels), and the bosonic elastic scattering rates (bottom 
panels), are depicted for single-beam (left) and crossed-beam (right) config- 
urations and a 6 Li- 23 Na mixture. The red-detuned beam is obtained with 
a Nd:YAG laser emitting at Ai = 1064 nm with a peak power Pi = 100 W. 
For jit > 400, the ratio between the blue- and red-detuned laser powers is 
kept at the constant value P2/P1 = 0.32. In the case of the single-beam 
configuration there is a significant decrease in the bosonic trapping fre- 
quency, while in the crossed-beam configuration a further gain through the 
increase of the fcrmionic trapping frequency is also evident. The large elas- 
tic scattering rates available during the entire evaporation process allows 
for a fast dynamics of thermalization of the Bose gas. We assume equal 
waists for the beams, w\ — W2 = 8 yum. The initial laser power is enough 
to have 77 = 10 for a 23 Na cloud consisting of = 10 6 atoms transferred 
from a magneto-optical trap at temperature T = 586 fiK (crossed-beam) 
or T = 222 fiK (single-beam), thermalized with a 6 Li cloud with Nf = 10° 
atoms. 
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mixture 


H (Hz) 


12 (Hz) 


7i b (Hz) 


7 2 b (Hz) 


"Li-^Na 


0.94 


12.5 


0.8 


63 


6 Li- 87 Rb 


0.94 


31 


4.7 


306 


40 K- 23 Na 


4.7 


7.8 


0.8 


63 


40 K- 87 Rb 


4.7 


9.4 x 10 3 


4.7 


1.1 x 10 4 



Table 3: Estimate of the Rayleigh scattering rates per unit of the corre- 
sponding laser power for the fermion-boson mixtures considered in the text. 
The laser wavelengths are the same as in Table 2 and the beam waists are 
chosen as w\ = W2 — 8 /im. 

In the following we will estimate the peak Rayleigh scattering rate, i.e. the 
value corresponding to the peak laser intensity in the center of the trapping 
potential, per unit of power of the laser beams in the case of the 6 Li- 23 Na 
mixture used in [3U| . The time dependence of the Rayleigh scattering rate 
can then be inferred by properly scaling the curves for the laser powers 
shown in Figure 7. As expected, the strongest source of heating comes 
from the blue-detuned laser beam acting on the sodium atoms. This is of 
less concern because the blue-detuned laser beam is turned on only during 
the latest stage of evaporation, and then raised to a fraction of the red- 
detuned laser beam. The residual Rayleigh scattering rate is estimated to 
be around 1.2 x 10 -2 Hz, much smaller than the estimated elastic collisional 
rate, and corresponds to lifetimes well in excess of 10 s. A similar situation 
occurs also for the 40 K- 23 Na mixture. As shown in Table 3, due to the 
proximity of the atomic transitions for 6 Li or 40 K and 87 Rb instead, the 
Rayleigh scattering rates due to the deconfming beam seems prohibitive. 
However, the large Rayleigh scattering rate per unit of power for the 40 K- 
87 Rb mixture is less frightening since this mixture has also a critical power 
ratio smaller by one order of magnitude with respect to the other ones, if 
equal waists are assumed. 

The lifetime estimated for 6 Li- 23 Na is long enough to perform various 
experiments aimed at evidencing superfluidity features, such as mechanical 
stirring or generation of collective excitations. 

5 Sympathetic cooling and possible evidences 
for a fermionic superfluid phase 

The study of the dynamics of evaporative cooling is a prerequisite to discuss 
sympathetic cooling of fermions through their interactions with the Bose 
gas. Heat exchange between two ensembles is perhaps the most widespread 
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thermodynamic process occurring in nature. Thus, as a general method to 
refrigerate an atomic or molecular ensemble whose direct cooling is diffi- 
cult to achieve, one can spatially and temporally overlap it with a colder 
ensemble, the so-called sympathetic cooling. This process has been first 
demonstrated at the microscopic level for trapped ions [71], then for neu- 
tral atoms and molecules via use of cryogenically cooled helium gas [75] , 
More recently, with the advent of ultracold atomic physics, sympathetic 
cooling in the microkelvin and nanokelvin ranges has been achieved for 
bosons in different internal states [7H|, in different isotopes [77], and for 
different species [75] , 

The elastic scattering properties among fermions and bosons at very 
low temperatures, which strongly influence the efficacy of sympathetic cool- 
ing, are starting to be collected for various atomic mixtures [30] E, ■ One 
potential problem common to all the mixtures is the diminished heat ex- 
change capability of bosons when approaching condensation. In particular, 
fermions can be considered as impurities in the boson cloud and below T c 
the Bose gas has a condensed fraction which is expected to be superfluid. 
Based on the Landau criterion for the critical velocity in a superfluid, we 
do expect Fermi-Bose collisions to be suppressed when the Fermi velocity 
becomes smaller than the sound velocity of the Bose gas |26|. This has 
been also experimentally demonstrated [7J5] • The Fermi velocity, defined 
by ZcbTf = can be expressed in terms of the trapping parameters 

as 

VF=1 . 91 ^y /2 ^/ 6 . (15) 

On the other hand, the sound velocity v s for the Bose gas can be written in 
terms of the chemical potential as /i = misleading, in the Thomas-Fermi 
limit, to 

where ab is the boson s-wave scattering length. The vp/v s ratio is, there- 
fore, 

/ h \ 1 / 10 /™ \ V2 /, , \ 1/2 w i/6 

vf _ v n \ ( "lb \ / Wf \ iv f 



— = 1.57 —\ — — (17) 

The following remarks are in order: a) the velocity ratio i|17|) scales 
with the square root of the trapping frequencies, so that it helps to have 
Wf/c^b >> 1 to avoid suppression of Fermi-Bose scattering; 6) the velocity 
ratio also depends upon the number of fermions and bosons. If the former 
are kept constant in the trap and the latter undergo evaporative cooling, 
the velocity ratio is increased; c) the ratio vp/v s is already large for con- 
ventional single-color optical dipole traps. For instance, if c^b = 2-7T x 10 4 
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s _1 and Nf — N\>, for a 6 Li- 23 Na mixture we have Vp/v s ~ 7.2. In this 
case, the loss of cooling efficiency becomes relevant for fermion velocities 
v < vp/7.2 corresponding to T/T-p < 2x 10~ 2 . It seems therefore that this 
mixture will hardly enter into the regime where superfluid suppression of 
impurity scattering is significant. 

A more stringent limitation to sympathetic cooling is instead set by 
the classical heat exchange between bosons and fermions. In Figure 8 we 
show the dependence of the ratio between the heat capacities for the Bose 
and the Fermi species, Cb/Cf , in bichromatic traps with single and crossed- 
beam configurations. This gives us a picture of the efficiency of sympathetic 
cooling since this process breaks down when the ratio Cb/Cf becomes of 
the order of unity. The heat capacities have been evaluated numerically 
as described in Appendix A taking into account the time dependence of 
the various parameters, in particular the diminishing number of bosons 
during forced evaporation. The effect of many-body interactions on the 
heat capacity of bosons, evaluated in |SU], is next-to-leading with respect to 
the dependence upon the number of bosons. We see that in the bichromatic 
traps considered here the heat capacity of the bosons is always larger than 
that of the fermions by an order of magnitude. 

In Figure 8 we also show the temperature ratios T/Tp and T/T c . In 
both single- and double-beam cases, it is possible to reach T/Tp < 1CP 1 
while for the Bose gas T/T c > 5 x 10 _1 , with some slight advantages 
in the crossed-beam configuration. This means that, with respect to a 
monochromatic optical dipole trap, a more substantial Bose thermal cloud 
can be sustained while the Fermi gas is in a deeper degenerate regime. In 
the various experiments reaching the Fermi degenerate regime, the estimate 
of the temperature is usually obtained by fitting the surviving normal Bose 
component. As the temperature is lowered, the thermal component shrinks 
in amplitude and size, and the fit to assess its temperature is less accurate. 
This effect is mitigated in our bichromatic traps thus allowing for a more 
precise thermometry. 

Simple estimates for the minimum degeneracy parameter T* /Tp can be 
obtained by extending the qualitative discussion reported in the concluding 
part of Section 2. Below T c the heat capacity of an ideal Bose gas can be 
written as |H3 



while for T/Tp < 0.5 with a good approximation the heat capacity of a 
Fermi gas is linearly dependent on the temperature (see also Figure 9) [S^j 




(18) 



Cf ~7T 2 fc B Af— . 



T 



(19) 
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Figure 8: Efficiency for sympathetic cooling of Fcrmi-Bose mixtures in 
bichromatic optical dipole traps. In the top panel the ratio between the 
heat capacities of bosons and fermions Cb / Cf is shown versus time for the 
evaporative cooling strategies chosen in Figure 7 for single-beam (dashed) 
and crossed-beam (solid) configurations. The time evolution of the tem- 
perature ratios T/Tp and T/T c is shown for the single-beam (center panel) 
and crossed-beam (bottom panel) configurations. 
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By taking the ratio of Cb and Cf and observing that according to Eqs. 

and J3J| T c can be expressed in terms of Tp, we obtain the degeneracy 
parameter T/Xp as 



In a conservative scenario we can assume that sympathetic cooling stops 
when Cb — Cf. In this case, for ui{/u)b ~ 1 (as in the case of the 6 Li- 7 Li 
mixture) we get a minimum value of the degeneracy parameter at the end 
of the cooling as T*/Tp ~ 0.35. If cooling is still possible for Cb < Cf 
more optimistic estimates can be given. For instance if cooling stops when 
Cb/Cf ~ 0.1, then T*/T F ~ 0.11. By using a larger u f /uj h ratio the T* /T F 
ratio decreases according to (20). It is interesting to note that the depen- 
dence is quite sensitive to the frequency ratio, and even minor deviations of 
this ratio from unity result in an observable effect. For instance, in a mag- 
netic trap with a 6 Li- 23 Na mixture we have Wf/cJb — (mb/nii) 1 ^ 2 ~ 1.96 
and therefore we estimate, for Cb/Cf = 0.1, a value T* /Tp = 0.04, very 
close to the minimum value T/Tf = 0.05 recently obtained by the MIT 
group for a complete evaporation of the Bose component |35) . 

Superfluidity of the Fermi gas is expected below the critical temperature 
for the onset of atomic Cooper pairs |23 ESI 



where fcp is the Fermi wavevector such that Ep = h 2 k F /2m{, and at is the 
fermion-fermion elastic scattering length. Besides leaving freedom to ap- 
ply arbitrary homogeneous magnetic fields to enhance the scattering length 
through tuning to a Feshbach resonance [441 1451 137] . as demonstrated ex- 
perimentally for Fermi gases in |561 157j , our bichromatic configuration al- 
lows also for an independent increase of kp due to the higher achievable 
densities. The resulting Tbcs/^f are within the accessible range which 
corresponds, as seen in Figure 8, to T/Tp >3x 10~ 2 . The use of opti- 
cal trapping also leads to a large absolute value of the Fermi temperature 
which can be otherwise obtained by magnetic trapping only with partic- 
ular geometries maximizing the field gradients |83|. Finally, the presence 
of the Bose gas could allow for enhancements of the BCS pairing tempera- 
ture since bosons can mediate phonon-exchange between fermions in a way 
analogous to ordinary superconductors |841 185) . 

In general, with the technique discussed above new mixtures consisting 
of a normal Bose gas (or a Bose condensed gas coexisting with a large Bose 
thermal fraction) and a degenerate Fermi gas are viable. This contrasts 
the only situation known so far of a 3 He- 4 He mixture where degeneracy 




(20) 



_2_ -7r/2fc F |a f | T 

3e 1 



(21) 
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is reached earlier for the bosonic species. One of the advantages of using 
such an anomalous mixture is the possibility to have a well controllable 
background - a normal Bose gas - superimposed to the Fermi gas. This 
considerably simplifies the possible signatures of a superfluid phase tran- 
sition in a Fermi gas. For instance it should possible to look at a bulge 
in the density distribution as predicted in since this is obtained ad- 
mist a smooth, low density and well controllable thermal cloud instead of 
a higher density and peaked condensate [SHj- The presence of a superfluid 
state could be evidenced also by using the same blue-detuned beam as a 
mechanical stirrer for the fermion cloud. In this case one should look at a 
finite threshold for the onset of a dissipative motion or a drag force [17U18) . 
The presence of the bosonic thermal cloud (or a condensate component) 
gives rise to heating for all stirring velocities (or above a critical velocity 
depending upon the condensate density) of the laser beam [HZj • However 
the contribution of the thermal cloud to the heating is very low, and much 
smaller than the Rayleigh scattering in the relevant range of stirring veloci- 
ties, due to the low density. In order to better discriminate heating coming 
from the Bose component, one could also take advantage of the proposed 
manipulation of an ultracold cloud with Raman beams |881 to create a 
directional critical velocity for the Fermi component. Another advantage 
of a Fermi-Bose mixture with the latter component in the normal state 
is the possibility to perform experiments on scattering from microscopic 
impurities (such as the one described in J75| for a Bose condensate) in a 
much simpler way than using two isotopes. 

The presence of the Bose cooler in the non-condensed phase allows one 
to re-examine also species for which Bose condensation has been proven dif- 
ficult to achieve, like cesium (see however the recent achievement of BEC 
in Cs with purely optical means |55j). For instance, one can reconsider the 
use of 133 Cs |5U] which, due to its large mass and small recoil temperature, 
can be efficiently cooled to very low temperature in a magneto-optical trap. 
This would ensure robust initial conditions, in terms of initial number of 
atoms and initial temperature, to start an efficient evaporative cooling in 
the optical dipole trap. Even if sympathetic cooling is less efficient due 
to the large mass ratio with the Fermi species [OT], the same feature is an 
advantage in terms of ratio between the trapping frequencies even at zero 
blue-detuned beam intensity or in a purely magnetic trap. Recently, sym- 
pathetic cooling of 7 Li through 133 Cs has been demonstrated in a far-off 
resonance optical trap A CO2 laser was used to create the optical 

potential obtaining, due to the static atomic polarizabilities, a larger en- 
ergy depth for 133 Cs. As a consequence the 7 Li went under simultaneous 
sympathetic and evaporative cooling. The problem can be circumvented 
by using a Nd:YAG laser as the primary, red-detuned trapping laser, and 
a deconfining beam with wavelength in between the two atomic transition 
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wavelengths. Furthermore, by relaxing the requirement for using an al- 
kali species (using for instance ytterbium cooled in an optical dipole trap, 
see |93l IM] , and recently brought into the degenerate regime for a bosonic 
isotope various favourable possibilities can be envisaged. 

6 Conclusions 

A novel path to reach deep Fermi degeneracy through sympathetic cooling 
with a Bose gas undergoing evaporative cooling in an optical dipole trap 
has been proposed. The key feature is that the trapping frequencies, deter- 
mining the degree of degeneracy of the dilute gases, are made different by 
properly using a second, deconfming beam. Both single- and crossed-beam 
configurations have been studied for the four Fermi-Bose mixtures available 
with stable alkali fcrmions. In all substantial increase of the Fermi 

energy is expected maintaining at the same time the critical temperature 
for Bose condensation constant or slightly decreased with respect to a single 
color optical dipole trap. This decrease of the critical temperature for Bose 
condensation is beneficial for maintaining a precision thermometry in the 
deep degenerate regime for the Fermi gas, and it makes also possible the 
use of species which can hardly reach Bose condensation. With respect to 
the routes for fermion cooling currently under active investigation, our pro- 
posal does not necessarily require enhancement of elastic scattering lengths 
through Feshbach resonances, a mechanism often associated to the increase 
of inelastic processes rates, in competition with hydrodynamic behaviour 
of the atomic clouds or formation of ultracold molecules [331 ESI CUD OH]- 
already present even in nondegenerate conditions With respect to 

the pure Fermi mixtures, it allows instead for a more precise determina- 
tion of the temperature and, in particular, of the BCS-like phase transition 
temperature below which onset of superfluidity is expected • Two tech- 
nical issues still to be addressed in detail are the degree of control of the 
laser beams to achieve a common focus, and the intensity ratio stability 
between confining and deconfming lasers. In a more pessimistic scenario, 
the proposed technique will allow to study the heating mechanisms pre- 
venting further cooling of fermions, in a way similar to the one proposed in 
|97|. or non-equilibrium phenomena in the ultracold regime analogous to 
those already explored for the Bose condensed gases 98 , or proposed as 
an interesting alternative to quasi-equilibrium sympathetic cooling |99| . 
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A Ideal Bose and Fermi systems with a finite 
number of particles in a harmonic poten- 
tial 

In this section we review some elementary thermodynamic properties of 
noninteracting quantum systems consisting of a finite number of bosonic 
or fermionic particles trapped into an external potential. More specifically, 
we suppose that the particles of mass m are confined by the anisotropic 
harmonic potential 

V(x, y, z) = l -ra (u%x 2 + u^y 1 + u> 2 z z 2 ) , (22) 
so that the corresponding single-particle energy levels are given by 

En x ,n v ,n„ = f^x (n x + ~J + tioJy (n y + + fuv z (n z + , (23) 

with n x , n y , n z = 0,1,2,.. .. 

In the grand-canonical ensemble, the average number of particles and 
the energy of the system are, respectively 
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n x =0n H =0n,=0 exp I ^ \ 
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where T is the temperature of the reservoir, /i the chemical potential, and 
the upper (lower) sign holds for fermions (bosons). Usually, one knows the 
temperature T of the reservoir and the number N of particles of the system 
so that the chemical potential must be evaluated as a function of N and 
T. The value of /i(N,T) is determined by solving the nonlinear Eq. I|24() 
numerically, e.g. by truncating the series in the r.h.s. at a proper order 
and controlling the convergence error |1U0| . Once fi(N,T) is known, the 
energy E(/j,(N,T),T) is evaluated in an analogous way by truncating the 
series in the r.h.s. of 125|) . 
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The computation of the series in 124fl or l|25ll requires the evaluation of 
triply nested sums which may be very time consuming. A more favourable 
situation, with double or single series to be computed, is obtained in the 
presence of symmetries of the external potential. For instance, in the case 
of experimental interest in which lu x = lu v = u> xy , Eq. I|24|) becomes 

oo oo 1 

" = EE r£.~-»ur (26) 

i».»=0n.=0 exp { — ^ — J ± 1 

where the factor n xy + 1 represents the degeneracy of the level at energy 
= tw xy (n xy + 1) + huj z (n z + i) with n xy , n z = 0, 1, 2, . . .. 
The heat capacity of the system at fixed number of particles is 



C(N,T) 



dE(N, T) 
dT 

dE dE dn 
df + 'dildf 
dE dE dN (dN 

7yr~~dJI Ijr \ d^ 



(27) 



where we used the constraint AN /AT = dN/dT + (dN/d(i)(dn/dT) = 0. 
The derivatives dN/dT, dN/dfi, dE/dT, and dE/d[i are obtained by JUJ 
and (|25|) and can be evaluated numerically as explained above. As an 
example, in Figure 9 we show the behavior of C(N, T) as a function of T 
for different values of the number of particles N. While in the fermion case 
C(N, T) is almost independent of N, for bosons the size effect is striking 
around the critical temperature. 

Finally, we derive the values of the Fermi and Bose degeneracy temper- 
atures reported in (JIJ and @- In the case of fermions, for any value of 
N and T we have n(N, T) < fi(N, T — 0) and the Fermi energy is defined 
as E F (N) = n(N,T = 0). For k B T > Hu, a condition well fulfilled in the 
experiments, the Fermi energy can be evaluated with the approximation 

N = £ 1 



±R3 F 3 

-fi^ (28) 



huj x fiWyhjj z Qtfiuj^ ' 

where tot = (uJxLUyLUz) 1 ^ 3 . Therefore, fcs?F = Ep = (6iV) 1 / 3 ?ia;f . 

In the case of bosons, we separate the occupation of the ground state 
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Figure 9: Specific heats of a system of N noninteracting bosons (left 
panel) and fermions (right panel) in a harmonic trap as a function of the 
normalized temperatures T/T c and T/T-p- We have chosen lo x — uj y — 
uj z /^/2 as in a crossed-beam optical dipole trap. 



from that of the excited ones: 
N = Nq + N e 

g k B T — 1 "x,n B ."z p k^T _ 1 

^0,0,0 

where fi(N, T) < E = E ,o, . When T -> 0, we have n -> £ and 7V 

The critical temperature T c is defined by the condition N e (fx = E$,T C ) — 

N. For k B T > ftw, we have 



poo 

N e (fi — E , T c ) ~ / dnsdnydriz 
Jo 



fc R T r \ 3 Z" 00 . , , 1 



dirdydz 



fiw h J J e x+y+z _ i 

C(3), (30) 



k-p.T, x " 



1/3 

where = (u> x oj y oj z ) and £ is the Riemann zeta function. Therefore, 
k B T c = (7V/C(3)) 1/: W 
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